Interplay of vacuum-mediated inter- and intraatomic couplings in a pair of atoms 
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The resonance fluorescence emitted by a system of two dipole-dipole interacting nearby four-level 
atoms inJ = l/2^J = l/2 configuration is studied. This setup is the simplest realistic model 
system which provides a complete description of the (interatomic) dipole-dipole interaction for 
arbitrary orientation of the interatomic distance vector, and at the same time allows for intraatomic 
spontaneously generated coherences. We discuss different methods to analyze the contribution of 
the various vacuum-induced coupling constants to the total resonance fluorescence spectrum. These 
allow us to find a dressed state interpretation of the contribution of the different interatomic dipole- 
dipole couplings to the total spectrum. We further study the role of the spontaneously generated 
coherences, and identify two different contributions to the single-particle vacuum-induced couplings. 
We show that they have a noticeable impact on the total resonance fluorescence spectrum down 
to small interatomic distances, even though the dipole-dipole couplings constants then are much 
larger in magnitude than the the single-particle coupling constants. Interestingly, we find that the 
interatomic couplings can induce an effect of the intraatomic spontaneously generated coherences 
on the observed spectra which is not present in single-atom systems. 

PACS numbers: 42.50.Nn, 42.50.Lc, 42.50.Gy 



a,: 

-I— > 
£ 

cr 



> 
o 
o 

CO 

7— I 

CN 

o> 
o 



I. INTRODUCTION 



Coherence and interference effects form the basis of 
many quantum mechanical phenomena [J ■ Its appli- 
cations have been revolutionized by the invention of the 
laser as a coherent source of light. But somewhat surpris- 
ing, under certain conditions it is also possible for coher- 
ences to be created in the interaction with the vacuum. 
An archetype system in which these spontaneously gener- 
ated coherences (SGC) have been predicted is the three- 
level V-type system shown in Fig. W[A-C) @, i, H|. The 
relevant physical processes can be understood intuitively 
as follows. Suppose the atom is initially in state |ei). The 
atom can be de-excited to the ground state \g) by the 
emission of a photon into the vacuum. This photon can 
leave the system, giving rise to spontaneous emission. Al- 
ternatively, it can be reabsorbed on the same transition, 
which leads to the Lamb shift. Finally, it could be reab- 
sorbed on the second transition, with final state \e2)- The 
latter process leads to the creation of SGC between the 
two excited states, and it has been shown in many theo- 
retical works that such SGC could give rise to fascinating 

applications 0, i, 0, 1 ! EI ED, HIM EI El, EJ, El El • 

The interpretation is facilitated by a quantized treat- 
ment of the light fields [l9| . In contrast, so far there 
are no conclusive observations of SGC in atomic sys- 
tems [2(| . The reason for this is that the process leading 
to SGC only occurs if both the two involved transitions 
are near-degenerate, and if their dipole moments are non- 
orthogonal. Speaking pictorially, the photon then cannot 
distinguish between the two transitions. Unfortunately, 
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these conditions usually are not fulfilled in real atoms. 
A proof-of-principle experiment verifying the presence of 
SGC could however be achieved in quantum dots [2lj . 

Recently it was found that there is a variant of such 
vacuum-mediated couplings with slightly relaxed condi- 
tions [H, [H, 01 ■ In this case, the two transitions in- 
volved do not share a common state, as it is the case 
e.g. with \g) in the V-type scheme. The simplest ex- 
ample is the four- level J=l/2<-»J=l/2 setup, see 
Fig. [UJD-F). One might be tempted to conclude that 
vacuum-mediated interactions between the two transi- 
tions are not possible, because two ground states of the 
transitions arc orthogonal, (S^ 5^) = with S± = |3)(1| 
and S2 = |4)(2|. But it turns out that there are ob- 
servables which nevertheless are affected by SGC. An 
example is the resonance fluorescence spectrum, which 
depends on (S^~ (f)S£" (t + r)). This two-time correlation 
was shown to be non-zero in general, which can be under- 
stood from the fact that the atom may evolve between 
the two ground states |3) and |4) in the time delay r. 
It is therefore possible to observe effects of SGC in re- 
alistic atomic systems. In the following, we will denote 
couplings of this type as intra-atomic couplings or single 
particle vacuum couplings (SPVC), since two transitions 
within the same atom are coupled. 

Similar vacuum-mediated photon exchange processes 
can also take place between two transition dipoles be- 
longmgto different atomsjp, [H, H3, El, H Hfl, [III, [S3, 
[H, MM HE M, [H, [H, EOT Fig. QIG-I) illustrates this 
for two two-level systems. Here, one atom is de-excited 
and emits a virtual photon, and the second atom ab- 
sorbs this photon and is excited from the ground state 
to the excited state. Processes of this type are known as 
dipole-dipole couplings, and will be called inter-atomic 
couplings of two-particle vacuum couplings (TPVC) in 
the following. The TPVC crucially influence the dynam- 
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FIG. 1: Vacuum-mediated couplings between different transi- 
tion dipole moments. (A)-(C) illustrate intraatomic or single- 
particle vacuum-induced (SPVC) couplings between different 
transition dipole moments in a single atom. (D)-(F) show a 
generalized form of these SPVC in a four-level J = 1/2 <-> 
J = 1/2 level scheme. (G)-(I) depict interatomic or two- 
particle vacuum induced (TPVC) couplings between transi- 
tion dipole moments in different atoms. In all three cases, 
the system is initially excited on one of the involved transi- 
tions, see (A,D,G). In the second step, a virtual photon is 
emitted together with the de-excitation of the atom (B,E,H). 
Finally, the virtual photon is reabsorbed on the second tran- 
sition (C,F,I). 



ics of the combined system, as can be seen, for example, 
from the resonance fluorescence spectrum [34|, ■ Such 
an energy transfer process between two particles is only 
possible, if the distance r separating the two atoms is 
small compared to the respective transition wavelength, 
and if the two transitions are near-degenerate. In con- 
trast to the SPVC, however, in general for TPVC there 
is no restriction on the dip ole moments of the two tran- 
sitions 0, [H, M, HI, M, El- The TPVC between or- 
thogonal dipole moments only vanish in certain relative 
alignments of the interatomic distance vector and the 
involved dipole moments. Therefore, in general TPVC 
have to be described using complete angular momentum 
multiplctts in order to obtain correct results [38| . 

From the above discussion it is clear that TPVC and 
SPVC are closely connected, and this is also reflected in 
their similar theoretical description. Motivated by this 
and by the many applications that have been suggested 
for TPCV and SPVC individually, here we study the sim- 
plest realistic atomic system in which both TPVC and 
SPVC can occur. This system consists of two dipolc- 
dipole interacting four-level systems inJ=l/2<-»J = 
1/2 configuration, see Fig. [2l In this setup, each in- 
dividual atom is modelled by a complete set of angu- 
lar momentum states, and fulfills the conditions for the 



generalized SPVC. In addition, the four dipole- allowed 
transitions in each atom interact with the corresponding 
four transitions in the second atom, giving rise to TPVC 
both between parallel and between orthogonal transition 
dipole moments. Our general aim is to study the impact 
of the couplings of the system dynamics and on the op- 
tical properties. More specifically, we are interested in 
the role of the SPVC, and the dependence of the SPVC 
contribution on the interparticle couplings. As main ob- 
servable, we discuss the resonance fluorescence spectrum 
of the two-atom system. In order to understand the ef- 
fects of the various couplings, we analyze different inter- 
particle distance classes. In the small-distance limit, the 
TPVC coupling coefficients are much larger than all other 
relevant system parameters. In the large-distance case, 
the TPVC vanish, such that the SPVC are the dominant 
vacuum couplings. In an intermediate case, both TPVC 
and SPVC coupling parameters are of similar order. We 
interpret our results using two different methods. The 
first analysis is based on the eigenvalues of the matrix 
governing the system dynamics, which describe position 
and width of the different eigenstates of the system. Sec- 
ondly, by artificially switching individual couplings on or 
off in the analysis, the quantitative impact on the re- 
sults can be studied. A combination of both methods 
allows to understand the formation of multiply dressed 
states in the system in detail. Regarding the SPVC, we 
find that the contributions to the spectrum by the single 
particle vacuum-induced couplings survive even at very 
low interparticle distance, where the TPVC coupling co- 
efficients are much larger than the corresponding SPVC 
coefficients. Therefore, also in the two-particle case, the 
J=l/2<->J=l/2 configuration enables one to ob- 
serve spontaneously generated coherences in a realistic 
level scheme. We distinguish two different types of in- 
traatomic couplings, which contribute either directly in 
the equations of motion, or in the total expression for 
the resonance fluorescence. We find that in particular 
for small interatomic distances, the parts entering the ex- 
pression for the spectrum are the dominating intraatomic 
coupling contributions, as in the single-particle case. But 
in contrast to the single-particle case, also the SPVC con- 
tribution in the equations of motion can significantly con- 
tribute to the obtained fluorescence spectra in the two- 
atom case. We thus conclude that the TPVC can induce 
additional SPVC. 



The outline of the paper is as follows. In Sec. Ill Al we 
provide the theoretical background for our analysis. We 
start by deriving the equations of motion for our biatomic 
system. We then derive expressions for the resonance 
fluorescence spectrum as our main observable. Finally, 
we present the eigenstatc spectrum of our system as a 
tool to analyze the contribution of the different coupling 
constants. Sec. IIIII contains our results. We start by 
discussing the eigenstate spectrum in Sec. ITTTAl and then 
proceed to discuss the role of interatomic couplings in 
Sec. IIII Bl Finally, in Sec. IIII C| we analyze intraatomic 
couplings. Sec. IIVI concludes with a summary. 
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FIG. 2: The analyzed model system consisting of two nearby 
four-level atoms in J = 1/2 « J = 1/2 atoms. This 
configuration is the simplest realistic system in which both 
single-particle and two-particle vacuum-induced couplings oc- 
cur. In (a) the geometrical setup is shown. We assume that 
atom 1 is located at the origin and atom 2 at the point V2- 
(b) shows the inner structure of each of our two atoms in 
J = 1/2 — * J = 1/2 configuration. Note that the energy 
differences are not to scale. 



II. THEORETICAL ANALYSIS 

A. Model system 

We consider a system consisting of N = 2 atoms in 
J = | <-> J = i configuration. Each atom has L = 4 lev- 
els and D = 4 possible electric dipole allowed (El) tran- 
sitions with dipole moments d i (i £ {1,..,4}) as shown 
in Fig. [2] d 1 and d 2 couple to linearly polarized light, 
the so-called it transitions, whereas the a transitions d 3 
and d 4 couple to circularly polarized photons. We define 
the mean transition frequency as 



UJ 



1 D 



(i) 



where u>i denotes the transition frequency of the i-th tran- 
sition. Calculating the dipole matrix elements for our 
system via the Wigner-Eckart theorem [42j we find 

d 1 = (l^e? 1 3),, =-^=De z , d 2 = d 1 



d 3 = (2| Al d|3) AI =yl^e_, d 4 =(d 3 )*. 

Here, T> is the reduced dipole matrix element, e_ = 
(e x — ie y )/y2 the circular polarization vector, /i g {1,2}, 
and e x ,e y ,e z are the Cartesian unit basis vectors. Wc 
assume the system to be driven by a monochromatic laser 
beam propagating in the y direction, 



E L (t) =£ L e l(kLr - MLt) e- 



(2) 



£l is the field amplitude, the wave vector, lu^ the 
frequency, and we choose the polarization e = e z such 
that the driving field only couples to the 7r-transition 
dipoles d x and d 2 . 

The free evolution of our two atoms is governed by the 
Hamiltonian 

Vat = H £ E S i+ S f- + S i- S <+) > ( 3 ) 

fi— 1 i— 1 



with <S>/^_ the atomic excitation operator for transition i 
in atom [1, and is the corresponding de-excitation 
operator. In Schrodingcr's picture the interaction with 
the laser field is described by 

2 4 

U L = - fi E E + H.c.) . (4) 

/^— 1 i— 1 

Here, the position dependent Rabi frequency is defined 
as 



«i(r) 



dj • e £ L 



(5a) 
(5b) 



Since d x and d 2 arc antiparallcl, wc can define 

fi(r) = 0e ifei r = n x (r) = -0 2 (r) . (6) 

In a suitable interaction picture with Aj = lol — uii the 
laser detunings, the full Hamiltonian then reads 

2 2 

W = -^EE + (£U(r l >)Sf + + H.c. 

/j,— 1 i—1 

(7) 

Finally, the system dynamics can be described by the 
master equation 

d t Qat(t) = tt[H, Qat{t)] + C a g at (t) + C 7 g at (t) . (8) 
in 



Here. 



N D D 

£a~Q«t(t)=i E EE^f 

l*,"=l t=l j = l 



S?\S»_,Q at (t) (9) 



is a modification to the coherent part of the evolution 
arising from the dipole-dipole coupling (TPVC) of the 
two atoms, as can be seen from the restriction /1 7^ v of 
the summation. It leads to the formation of collective 
dressed states, similar to the well-known symmetric and 
anti-symmetric collective states in two interacting two- 
level atoms [l[ . It can be seen that couplings between all 
four transitions in one of the atoms to all four transitions 
in the second atom are considered. The incoherent part 
is given by 

N D D 

C r § at (t) = - J2 EE^fWteO 



/L4,I/=1 i—1 j — 1 



(10) 



This term essentially contains three types of contribu- 
tions. For fj, = v and i = j, i.e., absorbing and emit- 
ting transitions are identical, the term describes the usual 
spontaneous emission with the rates ji of the individual 
transitions i, 



7i 



(11) 
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Second, if the energy exchange occurs between different 
transitions in the same atom (/z = v, but i j), then 
the corresponding process is a SPVC governed by the 
coupling constant 



di\\d., 



(12) 



Here, the normalized scalar product of the two dipole 
moments accounts for the fact that SPVC only occur 
between non-orthogonal transition dipole moments. Fi- 
nally, the terms with /x ^ v and i ^ j describe TPVC 
with coupling constants T^. The incoherent and coher- 
ent TPVC coupling constants van be expressed as 



ij 



and 



K) T Im X (7vK 



(dj T Rex(r^)d 

with the help of the tensor 
{kf 



(13) 



(14) 



Xp„(r) 
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In the following we assume the energy differences 5 = 
5 a = 0. This corresponds to the case where no external 
magnetic field is applied. Then the laser detunings Ai 
and A 2 become equal and we denote them by A. Addi- 
tionally, all transition frequencies become equal u>i = luq. 



B. Resonance fluorescence spectrum 

The total resonance fluorescence spectrum is given by 
the Fourier transform of the two-time correlation func- 
tion of the electric field operators Q 



e -*5r( £-(-) ( H) t + T ~ )E (+) (B, t) > st dr . 

J 

(15) 

In the far field zone the positive frequency part can be 
calculated as Q 



E {+) (R,t) ■ 
1 

ineoRc 2 



E 



'free 
fl—1 l—l 



(16) 



The negative frequency part E^ > (R, t) can be found by 
Hermitian conjug ation of E^ + \R,t). Here, R = RR 
denotes the position of the photon detector and t = t — 
R/c the retarded time with c the speed of light. We have 
also introduced the cross product factor ^(d il R) which 
describes the polarization structure of the emitted light, 



C(d it R) = R x (R x di) . 



(17) 



In Eq. (|16p . the first term denotes the positive frequency 
part of the free field. If the point of observation lies 
outside the driving laser beam, it does not contribute to 
the correlation function. 

Since we assume a polarization sensitive detector, it 
can detect photons emitted by a and n transitions sepa- 
rately. For this, we choose the observation direction for 
the 7r light in the x-y plane as R^ = (1, 1, 0) T /y/2. Then, 
products of the form C{d t , R^)^ (d j; R^) with i e {1,2} 
and j £ {3,4} vanish. Thus we can observe the photons 
emitted by the tt transitions separately. 

Ignoring retardation effects [41| , the spectrum for the 
linearly and circularly polarized light can be expressed 



S w (w) 



2 2 

-a* £ ]T 



x (Sf + (t + r)5/_(t) ) st e iko ^ r »- r ^ dr (18) 



Here. 



l/(4^ 0J Rc 2 ) 2 C(d 1 ,B 7r )C(di,« 7 r)- 



Throughout our analysis, we keep the positions of the 
detectors fixed, such that is a constant prcfactor 
which we will neglect in our numerical calculations. 

Next, we decompose the transition operators in mean 
values and fluctuations, 



sr ± (t) = (sr ± (t)ui + ssf ± (t) , 



(19) 



where 1 = li g) I2 is the product of the two spaces be- 
longing to the two atoms. Inserting this in Eqs. Q18p it 
turns out that the mean value contributions lead to the 
coherent part of the spectra whereas the fluctuations de- 
termine the incoherent part Q. The coherent parts of 
the spectra evaluate to 



(20) 



where 2" c ^ h denotes the total coherent resonance fluores- 
cence intensity. 

Since the coherent part of the spectrum only consists 
of a delta-peaked contribution at the driving laser field 
frequency, throughout our analysis we will focus on the 
incoherent part of the resonance fluorescence spectrum. 
In order to calculate the incoherent spectra we now define 
functions containing the Fourier transform of two-time 
averages as 



(21) 



The (uj) can be evaluated using the quantum regres- 
sion theorem [43l. |44T| . Using these functions, the incoher- 
ent 7r spectrum can be written as 



— 1 — l 



(<D) 



(22) 
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C. Eigenstate analysis 

It is well-known that the interpretation of the dynam- 
ical properties of a given quantum optical system often 
is facilitated by the introduction of the system dressed 
states. These dressed states arc commonly defined as 
the eigenstates of the interaction picture Hamiltonian Q . 
For our purposes, wc follow a similar strategy, but include 
the incoherent parts of the master equation as well, since 
they are crucial for the SPVC and TDPC effects studied 
here. We start by writing the master equation Eq. ||5J) as 
a vector-matrix equation 



d t Qat(t) = M Qat{t) . 



(23) 



Here, Q a t(t) is a vector containing the 256 different den- 
sity matrix elements describing our system, and A4 is a 
256x256 matrix describing the dynamics of the system. 
Next, we apply a unity transformation hi 



M diag = UMU 1 
QD = UQat 



(24) 
(25) 



leading to a diagonal matrix M-diag and a transformed 
state vector qq. The solution of the transformed master 
equation is given by 



Qo{t) 



Qd{0) 



We express the complex eigenvalues of Aidiag 



as 



J 3 ' 



(26) 



(27) 



with j G {1, . . . , 256} and Xh v 3 ^ ^- Then, the compo- 
nents of the solution for the transformed density matrix 
vector can be written as 



[Qo{t)] 3 



'.Qd(0)], 



>i)t 



Qd(0)} 3 • (28) 



Similar to the dressed state analysis, we find that nega- 
tive real parts of the eigenvalues lead to an exponential 
decay of qd (t) , whereas the imaginary parts result in an 
energy shift. Therefore, the Vj can be used to under- 
stand the origin of the different lines in the resonance 
fluorescence spectrum. 
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FIG. 3: Energy shifts of the system dressed states, evaluated 
as the imaginary part Vj of the eigenvalues = \i °f the 
matrix M. governing the system dynamics. The interatomic 
distance vector has direction given by 6 = 7r/2 and (f> = 7r/4. 
(a) shows the case of large interatomic distance ri2 = lOAyr 
in dependence of the Rabi frequency f2. The four branches 
clearly show the splitting of the spectrum into a central fea- 
ture and Mollow sidebands due to the AC-stark splitting of 
the bare states. The horizontal dotted lines mark the values 
for fl chosen for calculations in the later sections, (b) depicts 
the energy shifts as a function of the interatomic distance for 
fixed Rabi frequency Q = 107-n-. The three different cases 
of large, intermediate and small distance can clearly be dis- 
tinguished. The horizontal dotted lines mark distances for 
which we will present resonance fluorescence spectra in later 
sections. 



III. RESULTS 

A. Eigenstate analysis 

We start by analyzing the position of the system's 
dressed states in order to identify different parameter 
ranges of interest. Fig.[3^a) shows the imaginary parts Vj 
of all eigenvalues of the matrix M.diag in dependence of 
the Rabi frequency and at a large interatomic distance 
Tfj,v = r\2 = lOA^. In this case, the two atoms essen- 
tially act independently. Here, the imaginary parts of the 
eigenvalues lie around 0,^1^/77,. and ±217/7^. It should 



be noted, however, that not all potential positions indi- 
cated by the eigenstate analysis do lead to a significant 
peak in the spectrum. In fact, in [23j] was shown that in 
the spectrum of a single such a four level atom only peaks 
at 0, ±57/7^ occur. This Mollow spectrum we also have 
to expect for two independent atoms. In Fig. [3jb), the 
positions of the eigenstates are shown against the inter- 
atomic distance r\i for fixed Rabi frequency Q = 107^. 
Three different distance regimes can be distinguished. 
In the large- distance case which is already reached at 
r i2 ^ 0.2A7T, the eigenvalue spectrum is dominated by 
the Mollow structure already shown in Fig. [21a). In 
the small-distance case r\2 < O^A^, different eigen- 
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FIG. 4: Linearly polarized resonance fluorescence spectrum 
for different interatomic distances ri2. In curve (a) we chose 
ri2 = 0.08A^, in (b) ri2 = O.IA^, in (c) ri2 = 0.2\ n and in (d) 
ri2 = 10A,r. In all curves the Rabi frequency is = Wjtt, and 
the laser field is applied on resonance A = 0. The interatomic 
distance vector is aligned with 8 — n/2 and <f) — 7r/4. 

value branches with positions depending on the inter- 
atomic distance are formed. These can be attributed to 
the TPVC-induced energy shifts, with a rj~ 2 3 -dependcncc 
of the interatomic coupling constants. Finally, in the 
intermediate distance case O-ObX^ < fi2 < O^A^, the 
eigenvalues overlap such that a clear interpretation be- 
comes difficult. 

Next, we will present results for the spectra in these 
three different regimes of interatomic distances, and in 
detail explain the structure of the obtained spectra based 
on the different coupling mechanisms. 



B. Interatomic couplings 

In the following we study the influence of the inter- 
atomic couplings on the resonance fluorescence spectra. 
Since these parameters are negligible for large ryi we 
mostly choose smaller distances as compared to the wave- 
length \n . We will see that additional peaks in the spec- 
tra occur due to these couplings between the two atoms as 
observed when calculating the eigenvalues of the matrix 
A4 in Sec. Ill CI Artificially turning on and off some cou- 
pling constants will help us to get a better understanding 
of the peak positions and thus the relevant physical cou- 
pling mechanisms. We assume our atoms to be located 
in the x-y plane which means 9 — tt/2 and we choose 
<b = 7r/4. Therefore, the couplings between one n and 
one a transition vanish and thus no couplings between 
one 7r and one a dipole can occur. 



1. Spectra for large and intermediate interatomic distances 
We choose the Rabi frequency as f2 = 107^. 
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FIG. 5: Resonance fluorescence spectrum at small interatomic 
distances. In (a) ri2 = 0.04A-n-, in (b) ri2 = 0.046A^ and in (c) 
ri2 = 0.06A^. All other parameters are chosen as in Fig. [4] 

In Fig. IU we show the linearly polarized spectrum ob- 
served from the (1,1,0) direction. As interatomic dis- 
tances we choose ri2 = 0.08A W in curve (a), 7*12 = 0.1A,r 
in (b), ri2 = 0.2A,r in (c), and the large distance r±2 = 
10X n in (d). Since the interatomic couplings vanish for 
large distances, in this regime the atoms behave like two 
independent particles. Then for the n spectrum we ob- 
tain one peak at the laser frequency where Co = lo—lol = 
and two sideband peaks at ±251 as in the Mollow spec- 
trum, see curve (d) in the upper subfigure. The curve 
for r 12 = 0.2A W is still similar to that for large distances. 
Only the amplitudes of the peaks differs from the Mol- 
low spectrum for large distances. Thus the interatomic 
couplings can be considered weak for a distance as small 
as 0.2A,r, which is in agreement with our findings from 
Fig. [3J If we further decrease ri2 the peaks at u = 
and uj = ±20 become lower and additional peaks arise 
at higher frequencies |w|, see curves (a) and (b). In Fig. [3] 
we marked the values r 12 = O.OSA^ and O.IA^ by dotted 
lines. This regime of interatomic distances gives rise to 
rather complicated dynamics, and several Vj occur very 
close to each other. This is also reflected in the spectra, 
which is characterized by overlapping peaks. The reason 
for this is that the interatomic coupling parameters are of 
the same order of magnitude as the laser Rabi frequency 
O in this intermediate distance regime. 



2. Spectra for small interatomic distances 

In Fig. [5l we show results for interatomic distances 
chosen as (a) ru = 0.04A T , (b) ri 2 = 0.046A7T, and (c) 
7"i2 = O.OGAtt. In this regime of small distances we can 
see four clear sideband peaks in the spectrum instead of 
two in the Mollow spectrum. Since these sideband peaks 
move farther away from the center Cj — when decreasing 
ri2, they can be associated to the TPVC which increase 
in magnitude with decreasing distance. For the linearly 
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FIG. 6: Here, we show the ir spectrum for the distance ri2 = 
0.04A^ while we artificially set some interatomic couplings 
equal to zero. In each curve of the upper subfigure except 
curve (d) only the couplings of one Q a are on. In curve (a) 
the parameters of Q2 are turned on, in (b) the ones of Q>,, 
in (c) the ones of C/3 and in (d) we plot the spectrum with 
all couplings on. In the lower subfigure we turn one group of 
interatomic couplings after the other on. In curve (a) only 
coupling parameters of Q2 are on, in (b) the couplings of Q2 
and Q 3 , in (c) the constants of Q 2 , Gz and G4 are on and 
curve (d) depicts the spectrum with all couplings on. All 
other parameters are chosen as in Fig. 2] 



polarized light we still obtain a peak at uj = which 
becomes higher and narrower with decreasing interatomic 
distances. These spectra again correspond to the results 
of Fig. [31 One can see that the accumulation of imaginary 
parts of eigenvalues Vj for intermediate ri2 splits up into 
clear branches when decreasing the interatomic distance. 



3. Interpretation of the spectra 

In this Section we proceed to interpret the origin of 
the different spectral features found in the previous sec- 
tion. In order to do so, we artificially switch off parts of 
the interatomic couplings, and observe the change in the 
corresponding spectra. Note that in the following discus- 



sion, the SPVC are always kept in the analysis. For this 
procedure we divide our interatomic coupling constants 
into five groups Q\, G2, Qs, Qi and Q5. These groups 
are defined as follows 

g 1 = {r^ J n^|ie{i ) 2}andje{3 ) 4}} 

U{r^ i n^|*6{3,4} and j G {1,2}} 
= { couplings between one n and one a dipole } 
02 = {r^,nf;Kje{3,4}and^j} 

= { couplings between two different a dipolcs } 
Q 3 = {T^,Q^\ /x + v and i,j G {1, 2} and i ± j} 
= { interatomic couplings between 

two different tt dipoles } 
g 4 = {rf/, ftf/l ft jt v and i,j G {3, 4} and i = j} 
= { interatomic couplings between 
two equal a dipolcs } 
g 5 = {r%,Dg\ [i^v and i,j G {1, 2} and 1 = j} 
= { interatomic couplings between 

two equal w dipolcs } . (29) 

Qi and Q2 also contain some intraatomic couplings. 
However, these are zero because of the orthogonality 
of the respective dipole moments. Additionally, if both 
atoms are in the x-y plane, which means 9 = n/2, all 
interatomic couplings of group Q\ vanish. Since our sys- 
tem contains only two atoms and the couplings do not 
change when exchanging \i and v we use as superscript of 
the interatomic couplings inter instead of \xv. Our clas- 
sification is done in such a way, that within one group 
the absolute values of both real and imaginary part of all 
pmter or Qm*e^ respectively, are equal. This still holds 
if we choose 9 ^ -|. For each group Q a we define one pair 
of coupling parameters (r a , Q a ) as representative for this 
group. These definitions are made as follows 



p jointer 

1 1 - 1 13 

p jointer 

1 2 — I l A 



1 3 — 1 1 



r 4 
r 5 



jointer 
1 33 
jointer 
1 11 



Qi =n rater 



'13 



2 =n™ ter 

3 = n 

= 



inter 
12 

inter 
33 

inter 



(30) 



The absolute values \T l £ ter \ and |f^f er | of all pairs 
of coupling constants belonging to the same group are 
equal. 

In Fig. [51 our results for the spectrum emitted by the 
7T transitions are shown. We choose 7-12 = CLCMA^ and 
the Rabi frequency Q = lO'jn- Since our atoms are lo- 
cated in the x-y plane, all couplings of Q\ vanish. The 
upper subfigure shows spectra obtained when only one 
group of couplings is on and the others are set equal to 
zero. In curve (a) all parameters of Q 2 , in (b) the ones 
of C/3, and in (c) the ones of Q4 are on. Turning only 
coupling constants of Q5 on gives us the same spectrum 
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FIG. 7: Incoherent resonance fluorescence spectrum emitted 
by the tv transitions for the case where the couplings of <7 2 
and Qz are considered and in addition we progressively switch 
on the parameters of C/ 4 . In curve (a) the coupling constants 
of C/4 are zero, in (b) they are multiplied by 0.1, in (c) by 0.3, 
in (d) by 0.6 and in (e) by one. All other parameters are as 
in Fig. 11 



as in curve (b). Thus wc conclude that the interatomic 
interaction between equal 7r dipoles and between two dif- 
ferent ones has similar impact on the dynamics. This is 
not surprising since d x and d 2 are antiparallel. 

All curves show peaks at u> = and Co = ±2£1 due to 
the driving laser field as in the spectrum without inter- 
atomic couplings, see Fig. |4j In the spectrum plotted in 
curve (a) wc can sec additional peaks at Co w d= | 1 = 
±286. 57^. These peaks arise due to an energy splitting 
of the atomic levels caused by the vacuum-mediated cou- 
pling of the atoms. Curve (b) shows additional maxima 
at approximately Co « ± | O3 1 = ±91.67^. Looking closely 
one can see that the peak frequencies are a little bit 
higher than the coupling parameter. The reason for this 
likely is the influence of the Rabi frequency on the posi- 
tion of these outer peaks. This also holds for other peaks 
analyzed below. In curve (c) peaks at Co ~ ±103. 27^ 
occur which is equal to ±|fi4|. 

Now we continue by turning on the coupling constants 
one after the other. In the lower subfigure of Fig. [()] the 
resulting spectra are depicted. Curve (a) is the same as 
in the upper subfigure (only the plot range is changed). 

Curve (b) shows the ir spectrum where the couplings 
of Q-2 and C/ 3 differ from zero. We can see four peaks in 
addition to the triplet around Co = 0. These new peaks 
are located at approximately Co ss ±|^2| = ±286. 57^ and 
Co w ±|Sl3| = ±91. 67^. Note that the maxima around 
±286. 57^ arc hardly visible in the graph, since other 
curves show peaks at the same positions. 

In order to obtain curve (c) we turn the couplings of 
Qi on in addition to Q2 and Q3. Then the peaks at about 
Co « ±|r23| vanish and new maxima around Co = ±2007^ 
arise. Since this is not equal to the value of any coupling 
parameter, this peak cannot be interpreted straightfor- 
ward as the peaks in curve (a) and (b). In order to find 
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FIG. 8: Splitting of the atomic energy levels due to the inter- 
atomic coupling parameters for a small distance ri2 = 0.04A,r. 
For this interatomic distance the values of the coupling con- 
stants are ISI2I = 286.57^-, \fls\ = |f2s| = 91. 67^ and |fi4| = 
103.27^. 



the origin of these peaks around Co = ±2007^ we turn 
the couplings of Qa progressively on, see Fig. [7j Here, 
progressively means that we multiply the respective cou- 
pling parameters by a factor p £ [0, 1]. For curve (a) the 
couplings of C/4 are set zero (p = 0), in (b) p = 0.1, in 
(c) p = 0.3, in (d) p = 0.6, in (c) p = 0.9, and in (f) 
the respective couplings arc on which means p = 1. We 
can observe that the peaks at about ±|03| = ±91.67„- 
split up into two peaks. This splitting becomes larger 
for increasing values of the factor p. For p = 1 where the 
couplings of C/4 are on, the respective peaks are located at 
Co w ±117^ = |0 3 | - |^4| and Co w 194.87^ = |0 3 | + l^l- 
Note that the maxima at about ±117^ cannot be dis- 
tinguished clearly from the ones caused by the Rabi fre- 
quency of the driving laser field at ±207^ and the cen- 
tered peak at u = 0. During this level splitting when 
slowly turning on the couplings of CJ 4 the peaks at about 
±286. 57^ = ±1^2! do not change. This means that the 
couplings between equal a dipoles (f2 4 ) and between dif- 
ferent a dipoles (Q2 ) enter the dynamics of our system 
independently. By contrast, the level splitting of Fig. [7] 
shows us that the couplings between ir dipoles and be- 
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FIG. 9: Impact of the intraatomic couplings on the reso- 
nance fluorescence spectrum emitted by the 7r transitions for 
intermediate distances. The parameters are ri2 = 0.09A^, 
SI = 67^, A = —147^-, 6 — tt/2 and (/>=?. In curve (i) the 
complete spectrum with all couplings is shown as it could be 
observed in an experiment. In (ii), all intraatomic couplings 
are artificially turned off. In (iii), all intraatomic couplings 
entering the equations of motion are kept, while those enter- 
ing the expression for the spectrum Eq. (|18[) are set to zero. 
In (iv), the intraatomic couplings in the expression for the 
spectrum are kept while those in the equations of motion are 
dropped. 



tween equal a dipoles influence each other. 

Additionally turning on the parameters of Q§ we ob- 
tain the spectrum with all interatomic couplings as al- 
ready plotted in the upper subfigure of Fig. [6l Here, the 
peak at ±1^2 1 = ±286. 57^ no longer exists but is re- 
placed by maxima at w w ±(|fi2| ± |^s|) = ±(286. 57^ ± 
91. 67^). Note that for our parameters | £^3 1 + \Vti\ = 
91.67* + 103.2 77r = 194.8 77r « |0 2 | - |fi 5 | = 286.5 77r - 
91. 67* = 195. 57* and therefore the peaks at these posi- 
tions cannot be distinguished. To confirm the splitting 
of the peaks at Co w 286. 57* in curve (c) we do the same 
procedure for £7 5 as done for . We turn all couplings of 
G$ progressively on while the parameters of Q2, Q3 and 
Q4 are always on. Here we found a similar peak splitting 
as for £?4. 

Finally, we conclude our findings in this section by 
summarizing the various splittings in a dressed state 
analysis shown in Fig. [51 



C. Intraatomic couplings 

1. Classification of intraatomic couplings 

We now turn to the intraatomic (SPVC) couplings. 
In our system we can distinguish between two different 
types of intraatomic couplings. The first type are the 
coupling constants entering in the equations of motion 
for our system, F^f and 0^ ■ (where f2£t- are always zero) 



with i j. The second type are contributions of second- 
order correlation functions to the spectrum in Eq. (|22p 
where fx = v. We are thus led to split the resonance 
fluorescence spectrum into four contributions, 



C7T 



(w) = Vi(G>) + P a (w) - V 3 (u>) - Pa(w) , (31) 



where 



2 2 



V 2 (Cu) 



V 3 (u>) 



V±{Co) 



fl—1 l — l 

2 2 



E E r r(-): 



fi^u—l i—1 
2 2 



E E 



2 1 



E E r r(-)' 



(32a) 
(32b) 

(32c) 

(32d) 



fi,u—l i,j=l 



We then denote V\ — V3 as the intraatomic part of the 
contributions to the spectrum, and Vi + V2 — V& as the 
interatomic contribution. It should be noted, however, 
that this separation is difficult in particular due to the 
term ±7 , 4 . This term contributes for \i 7^ v which indi- 
cates an interatomic contribution, but at the same time 
it only contributes for i ^ j, which is the characteris- 
tic for an intraatomic coupling. One consequence of this 
is that spectra which are calculated by artificially sup- 
pressing e.g., the intraatomic parts, are not necessarily 
positive definite, which makes a straightforward inter- 
pretation difficult. It is important to keep in mind that 
results obtained by artificially suppressing parts of the 
couplings do not correspond to physically observable sit- 
uations, but can only be used to interpret the obtained 
full spectra in certain situations, as explained in more 
detail below. 



2. Effect of the intraatomic couplings on the resonance 
fluorescence 

We start our discussion by noting that in principle, 
both the intraatomic couplings in the equations of mo- 
tion and in the expression for the spectrum can crucially 
modify the total observed spectrum. This is surprising, 
since these couplings in the equations of motion do not 
contribute to the steady state, and thus the resonance 
fluorescence spectrum, of the corresponding single atom 
system [23| . The relevance of both types of couplings can 
most readily be seen for parameter ranges in which the 
intraatomic coupling constants T^ 1 are of similar mag- 
nitude as the interatomic coupling constants 0^- and 
^ij (m 7^ v )- An example is shown in Fig. [5] The 
figure shows the total resonance fluorescence spectrum 
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FIG. 10: Impact of the intraatomic couplings on the reso- 
nance fluorescence spectrum emitted by the n transitions in 
the small distance case. The parameters are as in Fig. [5ja) . 
The three subfigures separately show the three spectral fea- 
tures visible in the right half of the spectrum in Fig. [5ja). 
Each subfigure shows four curves (i)-(iv). As in Fig. [9j these 
correspond to the complete spectrum (i), the spectrum with- 
out intraatomic couplings (ii), and the spectrum with either 
the intraatomic coupling entering the expression for the spec- 
trum artificially suppressed (iii) or those entering the equa- 
tions of motion switched off (iv). 



(i) emitted by the it transitions, as well as the corre- 
sponding spectra obtained by artificially switching off the 
intra-atomic couplings in the equations of motion (iii), in 
the spectrum Eq. (TPS)) (iv), or both (ii). It can be seen 
that all four curves differ considerably, and we conclude 
that both types of intra coupling are of relevance. It 
is, however, difficult to attribute certain features of the 
observable total spectrum to either of the two contribu- 
tions. The reason is that for the parameters in Fig. [9l 
the total spectrum cannot simply be decomposed into 
the spectrum without intraatomic couplings and the two 
corrections arising from the couplings in the equations of 
motion and in the expression for the spectrum. 

We now turn to smaller interatomic distances in order 
to study the role of the SPVC in the case of numeri- 
cally dominating TPVC. Since the interatomic coupling 
constants then are much larger than the intraatomic cou- 
pling constants, typically the intraatomic coupling con- 
stants r^ 1 entering the equations of motion only slightly 



influence the total spectrum. An example is shown in 
Fig. [TQl for parameters as in Fig. [5]Ja). The three sub- 
figures show the three spectral features in the ranges 
[-50,50], [180,205] and [370,390] of the right hand part 
of the spectrum visible in Fig. 03a). Again, the total 
spectrum, as well as the corresponding spectrum with 
one or both of the intraatomic couplings artificially sup- 
pressed are shown. It can be seen that despite the clearly 
numerically dominating interatomic coupling constants, 
the spectrum is strongly modified upon suppression of 
the intraatomic coupling. In contrast to the example 
in Fig. O however, a grouping into two pairs of curves 
is observed. The curve for the full spectrum (i) coin- 
cides with the spectrum (iii) with intraatomic couplings 
entering the equations of motion artificially suppressed. 
Similarly, the spectrum (ii) without intraatomic coupling 
coincides with the result (iv) obtained by suppressing the 
intraatomic couplings in the expression for the spectrum. 
From this, we conclude that for these parameters, the in- 
fluence of the intraatomic coupling clearly arises from the 
corresponding parts in the expression for the spectrum, 
whereas the direct couplings in the equations of motion 
only give rise to minor corrections. This is the situation 
also found in the corresponding single atom system. 



3. Intraatomic couplings induced by interatomic couplings 

We found in Fig. [5] that already the intraatomic cou- 
plings in the equations of motion alone can give rise 
to a significant modification of the resonance fluores- 
cence spectrum, in contrast to the corresponding single- 
atom case. To interpret this difference in close analogy 
to the single atom first step, we evaluated 

the full steady state density matrix for parameters as 
in Fig. [9] both for the case with all couplings on, and 
with intraatomic couplings in the equations of motion 
suppressed. It turns out that while the magnitude and 
phases of the entries in the two density matrices differ, 
no elements are zero in one of the cases and non-zero 
in the other case. We thus conclude that no fundamen- 
tally new coherences are created due to the intraatomic 
couplings in the equations of motion. Next, for a bet- 
ter comparison with the single atom case, we analyze the 
steady state density matrix for one of the two atoms in 
our system. For this, we calculate the stationary state of 
the total system and then trace out the second atom. As 
a first example, Fig. [TTTa) shows the population of state 
|1) of atom 1 against the detuning A. From this figure 
it can be seen that at large intcrparticle distance, the 
intraatomic couplings entering the equations of motion 
do not have any effect, as expected from the single-atom 
case. However, at small distances, the populations with 
and without these couplings differ considerably in a range 
of negative detunings. In this sense, the intraatomic cou- 
plings between the two particles induce an effect of the 
intraatomic couplings which could not be observed in a 
single-atom system. This mechanism of interatomic cou- 
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FIG. 11: Intraatomic couplings induced by interatomic cou- 
plings, (a) shows the steady state population of state |1) in 
atom 1 obtained by tracing over the second atom, (b) shows 
the imaginary part of the coherence between states |1) and 
|3) in atom 1. The different curves are as follows: (i) All 
couplings included, large-distance case corresponding to no 
interatomic couplings, (ii) Large distance case, without in- 
traatomic couplings entering the equations of motion, (iii) 
and (iv) show the corresponding results for small distance 
ri2 = 0.09A^. The parameters are as in Fig. [9] except for the 
variable detuning A. Note that in (b), the two curves (i) and 
(ii) are shown multiplied with a factor 1/2 for better visibility. 

plings inducing intraatomic couplings is responsible for 
the dependence of the resonance fluorescence spectra on 
the intraatomic couplings in the equation of motion found 
in Fig. [5J In Fig. [TTT b). we show corresponding results 
for the imaginary part of the coherence between states 
|1) and 1 3) of atom 1. We again find that the interatomic 
couplings induce an effect of the intraatomic couplings 
entering the equations of motion. 

A possible explanation for these induced couplings is 
that in the single-particle case, the SPVC entering the 
equations of motion do not contribute to the steady state 
of the density matrix since the two ground states |3) and 
|4) are orthogonal, as explained in the introduction. In 
contrast, in the two-particle case, the relevant states are 
collective eigenstates originating from the TPVC, which 
each consist of different bare atomic states. Between 
these collective eigenstates, a modified set of transitions 
with different dipole moments occur. In particular, due 
to bare state mixing, near-degenerate non-orthogonal 
transition pathways originating from a single collective 
cigenstatc may be created. An analogous mechanism 



was suggested as a way to induce spontaneously induced 
coherences in a single three-level system in A configura- 
tion 0, [45[ . The two bare state transitions from the up- 
per state |e) to the lower states \a) and \b) originally are 
assumed to have orthogonal transition dipole moments. 
If a resonant laser field is applied to transition \a) «-> |e), 
then the system can be described in the dressed state 
basis {|+), R, |6)}, where |±) = (|e) ± |a))/V2. In this 
new basis, the two transitions from the dressed upper 
states |±) to the lower state \b) are near degenerate and 
non-orthogonal, such that quantum interference can take 
place. 

We thus conclude that the impact of the intraatomic 
coupling constants remains important even at low in- 
teratomic distances for which the interatomic coupling 
constants are much larger than the corresponding in- 
traatomic ones. For a large range of parameters, and 
in particular at smaller distances, the dominating contri- 
bution of the intraatomic coupling constants arises from 
the parts entering the expression Eq. (| 18[) of the spec- 
trum. In contrast to the corresponding single atom case, 
in a certain parameter range, the contributions entering 
the equations of motion can also have a substantial in- 
fluence on the obtained spectra. Since this contribution 
only occurs for the case of two nearby atoms, we conclude 
that the influence of the intraatomic couplings is induced 
by the intraatomic couplings. 



IV. SUMMARY 

In summary, we have analyzed a system of two dipole- 
dipole interacting nearby four-level atoms in J = 1/2 <-> 
J = 1/2 configuration. This is the simplest model system 
which on the one hand provides a complete description 
of the dipole-dipolc interaction for arbitrary orientation 
of the interatomic distance vector, and on the other hand 
allows for spontaneously generated coherences in a realis- 
tic atomic level scheme. The complete description of the 
dipole-dipole interactions is achieved by considering full 
Zeeman manifolds for the ground and the excited state. 
The spontaneously generated coherences contribute since 
this level scheme has two near-degenerate dipole allowed 
transitions with (anti-)parallel dipole moments. How- 
ever, as these two transitions do not share a common 
state, the comprehensive theoretical results on sponta- 
neously generated coherences such as in the usual V-type 
or A-typc configuration cannot be applied directly. We 
discuss different methods to analyze the contribution of 
the various coupling constants to the total resonance flu- 
orescence spectrum. A first analysis is possible based 
on the eigenvalue spectrum of the matrix governing the 
system dynamics. Then, we artificially suppress certain 
groups of couplings, in order to reveal their significance 
for the total spectrum. We also gradually switch on selec- 
tive couplings by artificially multiplying the correspond- 
ing coupling coefficient with a parameter ranging from 
zero to unity. This allows us to find a dressed state inter- 
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pretation of the contribution of the different interatomic 
dipole-dipole couplings to the total spectrum. 

Regarding the intraatomic couplings, we identify two 
different types of contributions. The first is via the 
SPVC-induced coupling coefficients directly entering the 
equation of motion. The second contribution appears in 
the expression for the resonance fluorescence spectrum. 
In general, both contributions can substantially influence 
the total resonance spectrum, even though it is difficult 
to attribute specific spectral features to either of the two 
contributions. However, in particular for smaller inter- 
atomic distances, we find that the dominant contribution 
of the intraatomic coupling is the one appearing in the ex- 
pression for the spectrum. A simple interpretation is that 



for small interatomic distances, the interatomic coupling 
coefficients entering the equations of motion are much 
larger than the corresponding intraatomic ones, such that 
they typically have only a small contribution. Neverthe- 
less, we find that the intraatomic couplings have an ob- 
servable impact even at small interatomic distances, such 
that also the two-particle J = 1/2 «-» J = 1/2 system 
is an interesting candidate to observe these intraatomic 
couplings. Finally, we could show that the intraatomic 
couplings entering the equations of motion can have a 
significant effect on the observed spectra, in contrast to 
the single-particle case. We thus conclude that the inter- 
particle couplings can induce additional contributions of 
the intraatomic couplings. 
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